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This paper is devoted to the study of the process e^e~ J/i/'X(3940) in the framework of 
light cone formalism. In our calculation two hypotheses about the structure of X(3940) meson are 
considered: X(3940) is 3^So state and X(3940) is one of 2^P states. The former hypothesis leads 
to a good agreements with the cross section measured at the experiment. As to the latter one, it is 
proposed a mechanism that allows one to understand the suppression of 2P mesons production in 
hard processes. 
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I. INTRODUCTION 



There are a number of charmonium-like mesons discovered recently in different experiments (see review 0). Our 
aper is devoted to a charmonium-like meson A"(3940) discovered in inclusive process e~^e~ —f J/ ■0+ anything at Belle 
The distribution of masses recoiling against the reconstructed J/^jj in e~^e~ J/'i/'+anything events measured at 
this experiment is shown in Fig. 1. 
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Fig. 1: The distribution of masses recoiling against the reconstructed J/ip in e'^e~ J/?/)+anything events. 



Evidently the theory that pretends to the understanding of the mechanism of A'(3940) meson production in e+e"- 
annihilation must also describe well the processes e+e" J/tjjr]c, J/4'Vc> J/'4'XcO measured at Belle. Considering 
charmonium mesons as a nonrelativistic bound states of cc pair one tried to apply NRQCD to these processes 3, but 
the results of these calculation proved to be in contradiction with Belle and BaBar measurements 0, [I]. 

Another approach to the prediction of the cross sections of these processes is light cone formalism. The authors of 
paper Q calculated the process e~^e~ — + J/iprjc in the framework of this approach. Despite of the uncertainties the 
agreement with the experiments was rather good. Further progress in understanding of exclusive double charmonium 



production in e~^e annihilation was connected with papers and Q. In the former paper the process 



was calculated. The latter one predicts the cross section of the process e+e" —i- J/^XcO- The values of the cross 
sections obtained in these papers are in good agreement with experimental results. In addition to good agreement 
with the experiments the application of light cone formalism allows one to understand that the wave functions of 
charmonium mesons are too wide to describe double charmonium production in the framework of NRQCD. 
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In our paper we will consider X(3940) meson in the framework of conventional quark model and sup- 
pose that this meson is an excitation of charmonium meson already seen at Belle experiment. There are 
tree mesons r]c{l^ So),ri'^{'2^ So),Xco{^^Po) in Fig- 1- Actually the peak labeled by Xcoi^^Po) can arise from 
Xco(l'^-Po), Xci(l'^-F'i), Xc2(l'^^2) mesons, although Xco gives dominant contribution. So two scenarios must be consid- 
ered: X(3940) is Ti'^iS^So) state and X(3940) is one of Xco(2^-Po), Xci(2^^i), Xc2(2^^2) states. 

This paper is organized as follows. In section 2 the model for the light cone wave functions used in our calculation 
is considered. Section 3 is devoted to the calculation of the cross section of X(3940) meson production. In last section 
we discuss results obtained in our paper. 



II. THE LIGHT CONE WAVE FUNCTIONS OF ^So,^ Si AND ^Po THE MESONS. 

In this section the light cone wave functions of ^S'o^iSS and '^Pq mesons will be considered. The functions used in 
our calculation are defined as follows (for details see 
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Here the following designations were used: Mq — Mq^'{ij = Mq^), x = xi,X2 = 1 — a; are the fractions of meson 
momentum carried by c quark and c antiquark correspondingly, the factors Zp, Zt, Z.„n Z^ are defined as 



Zt - 
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where Cp — 4/3, ho — 25/3. The wave functions (pi = Vp, V±, Vr, Va, Pa, Pp, Ss are normalized as follows: dxcjji — 

1. The wave function Sy is normahzed as dx{xi — X2)Sv{x) = 1. The dependence of the light cone wave functions 
on the scale is very slow and it will not be considered in the full form in all functions used in our calculation. Only 
renormalization factors of the corresponding local currents Zp, Zt, Zy will be regarded. 

Unfortunately today there is no information about the light cone wave function obtained directly from QCD 
Lagrangian. So to proceed with numerical analysis we are forced to use some model for these wave functions. To find 
the leading twist wave functions Vl{x), Vt{x), Pa (a;), Sv{x) we will apply Brodsky-Huang-Lepage(BHL) procedure 
which allows one to connect the light cone wave functions of leading twist with the equal time wave function in the 
rest frame. The equal time wave functions of charmonium mesons will be taken from the solution of Schrodinger 
equation with BuchmuUer-Tye potential 10]. Having these wave functions in momentum space ip(k^) one can get the 
light cone wave functions of leading twist using the following rule [llj : 



d'^k±il^cix,'k±). 



where ^pc{x,'k±) can be obtained from V'(k^) after the substitution ^ 
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Here nic is the quark mass in the potential model. 




Fig. 2: The functions ^{x) for IP, 2P state mesons. 

It is worth noting that in paper the relations between the hght cone wave functions and equal time wave functions 
of charmonium mesons in the rest frame were derived. The procedure proposed in paper is similar to BHL with the 
difference: in formula (O one must make the substitution (Pk± d^k^^J^^+ml/ {AmcXiX2)- But this substitution 
was derived at leading order approximation in relative velocity of quark-antiquark motion inside the charmonium. At 
this approximation ^ 0{v^), AxiX2 ~ 1 + 0{v'^) and the substitution amounts to dPk±^ (Pk±{l + 0{v'^)). Thus 
at leading order approximation applied in 12] these two approaches coincide. 

Using equation (O the expression for the leading twist wave functions can be easily written in the following form 



(7) 



where the constants Ci are fixed from the normalization condition, (pf^ix) are the asymptotic forms of wave functions 
of leading twist ^, the functions ^i{x) are given by the formulas 
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for P-wave mesons. 

For the light cone wave functions of nonleading twist there is no relation similar to equation ([5]) and to calculate 
these functions the following model will be applied. It is known that nonleading twist wave functions can be written 
as a product [ll| 



Ux)^<i^'{x)~^,{x) 



(9) 



with unknown functions (x) . In our calculation we will suppose that the functions <I>i (x) equal to the corresponding 
functions $i(x) ([7]) of leading twist. Thus the model for the light cone wave functions of leading and nonleading twist 
is given by equations O-®. 



III. THE STUDY OF THE PROCESS e+e" ~+ J/*X(3940). 

First let us consider the hypothesis: X(3940) is one of Xco(2^Po), Xci (2^Pi), Xc2 (2^P2) mesons. At Belle 
X(3940) is seen to decay to DD*^ and not to DD. If this DD* is dominant decay mode than X(3940) is Xci(2^-Pi). 
Unfortunately the process e+e~ J / il'x'cii'^^ P^) has not been considered in the framework of light cone formalism 
yet but this formalism tells us that the cross section of the process e+e" — > ^/V'X^i(2^-Pi) is suppressed by the factor 
1/s in comparison with the cross section of the process e+e^ — > J/'(/'x^q(2^Po) y. So one can expect that the cross 
section of Xco production is greater than the cross section of x'ci production and if x'ci is seen at the experiment Xco 
meson must be seen also. But the decay mode of x'co meson is DD and since this decay mode is not seen at the 
experiment we reject this hypothesis. 



^ for the wave functions Vl,V^,Vt ,Va-, Pa, Pp , Ss f/""" ~ 2:1x2, for the wave function Sy 'p"'" ~ xiX2{xi — X2). The asymptotic forms 
of the other functions used in the calculation can be found in papers @, 0, HI 




Contrary to the NRQCD predictions the values of the cross sections of exited charmonium mesons production 
measured at Belle and BaBar experiments are rather large. For instance due to this the process e~^e~ J/'^j]'^ with 
excited j]'^ meson production is seen at the experiments. In connection with this fact the question arises: if the process 
e^e^ — > J/'i'XcO is seen at the experiment and X(3940) is not x'co meson why one does not see the production of 
excited charmonium state x'cO- of the possible answers to this question is presented in Fig. 2 where the functions 
<i>(a;) ([5]) for IP and 2P states are shown. From this plot one can see that contrary to IP state the function of 2P 
state ^(x) is oscillating function. In light cone formalism this wave function must be integrated with the hard part 
of the amplitude of the process, consequently one can expect considerable cancellation in the amplitude. We will not 
estimate the value of the cross section of x'^q meson production since the accuracy of our calculation for this process 
is very low. This can be understood as follows if we divide the amplitude of x'co meson production into two parts 
where functions ^2P is positive and where this functions is negative than the amplitude is the difference between two 
large close numbers. Evidently the error of such calculation is large. It should be noted here that if our conjecture 
is correct and the suppression of x'cO production really results from the oscillation in the wave function of this meson 
than one can expect that the production of x'cQiX'c1tX'c2 mesons is suppressed in any hard process. 

Now let us consider the hypothesis: X(3940) is excited ^ meson. Two diagrams that give contribution to the 
amplitude of the process under consideration are presented in Fig 3. The other two can be obtained from the depicted 
ones by charge conjugation. The amplitude of the process e"'"e~ — > J/ipP for final mesons with equal masses was 
first found in paper [6| . In our calculation the result of paper Q will be used where the case of different masses was 
considered. The cross section of the process involved can be expressed through the formfactor F^p defined as follows 

(y(pi,A),P(p2)|J^|0> - e^.pae>?p^i^„p, (10) 

The formfactor Fyp equals 
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where ~ (s — My — Mp), Pa, Pp, Vt, Vl, V±, Va are the light cone wave functions defined above, Mv, Mp are the 
mass of the vector and pseudoscalar mesons correspondingly, d{x, y), s{x), s{y) are defined as follows: 

d{x,y) = — == xiH h/iH , = (13) 
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Fig. 4: The functions $(a;) for IS, 2S, 3S state mesons. 

From formula (jl3p one sees that gluon propagator d{x,y) tends to infinity when the momenta of quark (ki) or 
antiquark (li) created by this gluon tend to zero. Evidently this property has nothing to do with real situation: 
d{x,y) AmI/qQ when x,?/ — > 0. Similar problem is seen in the expressions for quark propagators s{x),s{y). This 
problem appeared since the authors of paper 6] used the following expression for the quark momenta ki — (+, J- 
, — ) = {(IqX\, 0, Mq/xiqi)) ( and similarly for three other quarks). It is not difficult to understand that this expression 
can be applied if the energy of the quark is ^/s. If a; — > the energy of the quark is ~ q^xi <C -y/s and this expression 
becomes inapplicable. Nevertheless the calculation of the cross section of IS*, 25* charmonium mesons production can 
be done using equations (|10p - (|14p since these charmonium mesons can be considered as nonrelativistic objects. This 
means that the regions where expressions (1131) . (|14p become incorrect are suppressed. 

Contrary to the 15 and 25* charmonium mesons there is considerable contribution to the cross section 35* charmo- 
nium meson production from the regions where the expressions (1131) . (jl4l are incorrect. Thus these expressions must 
be modified. We will do this as follows. First we will disregard transverse motion of quark-antiquark pair even in the 
end point region x,y ^ since account of the transverse motion is higher order v effect. Next if the energies of all 
quark q^Xi, q^yi in Fig. 3 are greater than rric than the following expression for the quarks and gluon propagators will 
be used: 

d{x,y) = xiyi, 

s{x) = Xl, 

s{y) = Vi- 

If the energy qoXi of quark fci reaches the value rric, what means that this quark is at rest, the expressions of the quark 
momentum ki — {qoXi, 0, m^/xiqo) must be substituted by {rric, 0, rric) (and similarly for three other quarks). We will 
not write explicit expressions for gluon and quark propagators d{x,y), s{x), s{y) since they are rather cumbersome. 
The model that will be applied in our calculation is rude but it is sufficient to apply it for the estimation of the value 
of the cross section. It is worth to note here that in the calculation of quark propagators s(x),s{y) we distinguish 
running mass Mq{ii) inside the quark propagators from the pole mass rUc of the external quarks ki,k2,li,l2- Moreover 
in our calculation all scale dependent quantities will be taken at scale ~ s/4. 

There are two contributions to the formfactor F^jp factored in formula (jlip . The first contribution originates from 
the wave function of mesons at the origin and it is proportional to ~ fvfp- The second contribution regards internal 
motion of quark-antiquark pair inside mesons and it is proportional to Iq. The leading NRQCD approximation does 
not take into the account the contribution of the second type. 

In the numerical calculation one loop expression for with Aqcd — 200MeV will be used, Mq = 1.2GeV, 
rric — 1.4 ± 0.2. The functions ([8]) of IS, 2S, 3S meson states at /i ^ Mq are presented in Fig. 4. The constant 
fv — 0.4GeV is determined from the decay width T{J/ip e+e^). Supposing that X{3940) is 3^5*0 meson the 
constant fp can be estimated as fp ~ /y(3^5i). In turn the constant /v'(3"^5i) can be found from the decay width 
r('0(4O4O) e+e~). Thus we have fp ~ 0.17GeV. These numerical parameters lead to the value of the cross section 

cr(e+e~ ^ J/'(/'X(3940))~ 11 ±3 fb. (15) 

Experimental result for this cross section is cr x Br(X(3940) > 2charged particles) = (10.6 ± 2.5 ± 2.4)fbt25 and it is 
in good agreement with our prediction. Sure there are a number of uncertainties connected with different reasons. 
For instance, one does not know the size of QCD radiative corrections and the size of 1/s corrections. In addition to 
these uncertainties there is an error due to the model of the light cone wave functions ([7]) used in our paper. 

To calculate equal time charmonium wave functions nonrelativistic potential model was applied. Our calculation 
shows that there is considerable contribution to the cross sections from the kinematical region where the motion of 
quark-antiquark pair inside charmonium cannot be regarded as nonrelativistic. To make our prediction more reliable 
in addition to nonrelativistic potential model we have calculated the cross section in the framework of the model 
where charmonium mesons are treated as a relativistic quark-antiquark bound state 13]. In the framework of this 
model the value of the cross section is cr ~ 10 ± 3 fb. 
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IV. DISCUSSION. 



In our paper light cone formalism was applied to the process e+e^ J/^X(3940). We considered two hypotheses 
about the structure of X(3940) meson: X(3940) is one of Xc0iXci)Xc2 mesons and X(3940) is r?" meson. Based on 
experimental data we rejected the first hypothesis and proposed the conjecture why Xc mesons are not seen at the 
experiment. If this conjecture is correct we predict the suppression of x'c mesons production in any hard process. As 
to the second hypothesis we calculated the cross section of X(3940) meson production if this meson is 77". The result 
of our calculation is in good agreement with the experiment. 

There are a number of uncertainties of our calculation connected with different reasons. For instance, one does 
not know the size of QCD radiative corrections and the size of 1/s corrections. But we believe that the main source 
of uncertainty is our model for the light cone wave functions. The main problem is that in our calculation we used 
nonrelativistic potential model. At the same time there is considerable contribution to the amplitude under consider- 
ation from region where the motion of quark- antiquark pair inside mesons cannot be considered as a nonrelativistic. 
Obviously application of nonrelativistic potential models to this region results in large error. Nevertheless we believe 
that if one was able to build correct description of chamonium in the relativistic region the value of the cross section 
(fTS)) would be of order of 10 fb. So the estimation of the cross section ([15]) using nonrelativistic potential model is 
rather good. 

Last statement is based on the following facts. First of all in addition to the nonrelativistic potential model we 
applied relativistic potential model [ll| to the calculation of the cross section of the process e+e" J/^X{39'iO). 
In the framework of this model we obtained the value a ^ 10 ± 3 fb. Moreover it is possible to estimate the cross 
section in a model independent way. As was noted above the amplitude of the process e~^e~ —> J/'^P is a product 
of NRQCD result for this amplitude and the factor that regards internal motion of quark antiquark pair inside the 
mesons. For the cross section this statement can be written in the following form 

a{e+e- -.J/^P)^fp\Io{P)f. (16) 

Obviously the wave functions of higher charmonium states become wider. Consequently the factor Io{P) that regards 
internal motion of quark antiquark pair inside the mesons is larger for higher charmonium states. So /o(3S') > /o(2S'). 
Using this relation and equation ^TE\i one can get lower bound for the cross section of 77" meson production 

a(e+e- ^ J/* Vc ) > ^(e+e- ^ J/^ if/-^ ~ 6 fb. (17) 

Another estimation of the cross section can be obtained if one supposes that /o(2S')//o(lS') ^ /o(3S')//o(2S'). This 
estimation gives cr 8 fb. 

As it was noted earlier the motion of quark antiquark pair near to the end point regions(a; '--^ 0, x ^ 1) is relativistic 
and cannot be calculated reliably in the framework of nonrelativistic potential models. Recently in paper [l3| it was 
proposed to solve this problem in the framework of NRQCD. We are not going to discuss this paper in detail. We 
would like to say only that the results of |13 were not applied in our paper since we don't think that it is possible to 
regard relativistic motion as it was done in |12l | . 

The interpretation of X(3940) as 3^ So state was proposed in paper [l3|. The problem with this interpretation 
consists in the fact that the mass of this meson obtained in the framework of potential approach is 4040 — 4060MeV ij. 
But as it was shown in our paper relativistic motion of quark antiquark pair is important for 3^ So meson and probably 
potential models do not regard this motion correctly. Moreover in paper 1^] it was shown that X meson with the 
mass 3940MeV agrees well with Regge trajectory with quantum numbers ^. 
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